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Abstract 

We study Laplace- type operators on hybrid manifolds, i.e. on con- 
figurations consisting of closed two-dimensional manifolds and one- 
dimensional segments. Such an operator can be constructed by using 
the Laplace-Beltrami operators on each component with some bound- 
ary conditions at the points of gluing. The large spectral parameter 
expansion of the trace of the second power of the resolvent is obtained. 
Some questions of the inverse spectral theory are adressed. 



1 Introduction 



The spectral theory on compact Riemannian manifolds has been studied for 
a long time and takes its roots in physical problems. A great number of 
important results has been obtained and this subject has a lot of ramifica- 
tions. The main objects of investigation in spectral theory are Laplace type 
operators on a compact manifold, constructed from a Riemannian metric. 
A related subject which is developing very actively is the spectral theory on 
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manifolds with singularities, see e.g. [HE]. Generalizing in another related 
direction, one can be interested in the spectral theory on metric graphs, or 
quantum graphs [9], or more general multistructures pQ. In the present 
paper we are studying some spectral properties of the Laplacian on spaces 
composed of one- and two-dimensional pieces. 

We consider the so-called "hybrid manifolds" , which are unions of smoth 
manifolds connected by line segments. More specifically, let us describe it 
from the topological point of view. Consider a set of m compact Riemannian 
manifolds Mi, M m and a set of n segments L±, L n . On each manifold 
Mi we fix some points qi s , s = 1, . . . , fa, fa > 0, i = 1, . . . , m. First of all we 
consider the disjoint union of all initial elements: M\U- ■ -UM m UL\U- ■ -UL n . 
Then we construct a one-to-one correspondence between the set of end points 
of all segments and the set of points qi s , s = 1, . . . ,fa,i = 1, . . . , M. The 
following natural condition on the number of elements must be satisfied 

m 

S^i = 2n. 

i=i 

Finally, according to this correspondence, we glue each end of each segment 
to the corresponding point on one of the manifolds, which results in a certain 
topological space. We assume it to be path connected which immediately 
implies that n ^ m— 1. The topological space obtained by gluing the initial 
manifolds and segments as described is called a hybrid manifold. If the 
manifolds are zero-dimensional, then we have a quantum graph. However, 
from now on, we will assume that our manifolds are 2-dimensional. We will 
also often relabel the gluing points as gi, . . . q^, where 

N = 2n. 

Our first task will be to define Laplace type operators on a hybrid space. 
This construction involves the Laplacians on each smooth part as well as 
boundary conditions at the points of gluing, and will be explained in Section 
3. The spectral properties of the operators obtained in this way can be 
studied using their resolvents or some function of it. In our approach we 
consider the trace of the squared resolvent. This is a standard procedure 
in spectral theory. Indeed, if we consider the usual (positive) Laplacian A 
on a Riemannian surface M, then it is known (see e.g. [TB]) that R 2 (z) = 
(A + z 2 )~ 2 is trace class. If G 2 (x, y; z) is the integral kernel of this operator, 
then we have the following asymptotic expansion as z — > oo : for all K ^ 0, 



<?(*■*;*) = E M ^ +1 > +o(.-<«*>), (i.i) 



k=0 



where the coefficients at are the local heat kernel invariants of M and can 
be expressed by some universal expressions containing the metric tensor and 
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its derivates. Furthermore, 



T ^ 2 (-) = E^3^+^- 2( " +2) )' a- 2 ) 



where the coefficients 



47rz 2fcH 
fc=0 



at = / a fc (x) 



are the global heat kernel invariants of M. In particular, the first terms of 
this expansion reveal important geometric and topological characteristics of 
the surface: 

TrR2 vol M X (M) 

where volM and x(M) are the volume and the Euler characteristic of M, 
respectively; this is closely related to the Weyl asymptotics for the eigen- 
values, see e.g. |4]. Note that traces of some functions of the Laplacian 
are also of use in the inverse spectral theory of quantum graphs, see e.g. 

[annua nz]. 

Our aim is to obtain similar results for a Laplace type operator on a 
hybrid manifold. In fact, due to the singular structure of the hybrid space, 
we won't get an asymptotic expansion in the usual way but in a certain 
extended form. More precisely, we show (subsection I3.3[) 



Theorem 1. Let R{z) denote the resolvent of a Laplace type operator on a 
hybrid manifold. Then Tri? 2 (z) has an expansion of the following form: for 
any k > 4 and a £ (0, 1) one has 

TttfM = Ei^Mj JU c,-(lnz 2 ) Q (J_\ 

H[) 4vrz2 + 4z3 z k + u { z k+a)> 

=4 

where volMj (resp. L) is the volume of the manifold Mj (resp. of the 
segment Lj) used to construct the hybrif manifold, and the coefficients Cj 
are rational functions. 

Each coefficient c,- is actually a rational function of nonpositive degree. 
Therefore, each Cj(ln z 2 ) has an aymptotic expansion in powers of 1/ In z 2 for 
large z. For a rather large class of boundary conditions (see the beginning 
of Section 5 for more details), we are able to compute the first terms of the 
expansion of c&. This allows us to address some quesions in inverse spectral 
theory : given the expansion of TrR 2 (z) in Theorem [H can we recover 
some topological and metric parameters of the hybrid manifold from the 
coefficients? Under some generic assumptions we prove the following result 
(subsection | 
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Theorem 2. Consider the expansion of the trace of the square of the resol- 
vent as in Theorem^ The knowledge ofTrR 2 determines: 

• whether this manifold is hybrid or smooth; 

• the sum of the volumes of all manifolds taking part in the construction; 

• the sum of the Euler characteristics of all manifolds; 

• the number of segments used in this hybrid manifold; 

• the sum of the lengths of these segments; 

• the Euler characteristic of the hybrid manifold. 

We mention that a rather similar result for quantum graphs was obtained 
in [12]. 

Finally, we will show in Theorem [22] that the expansion given in The- 
orem Q] contains some information about the boundary conditions used to 
construct our Laplace type operator. We note that there is a certain canon- 
ical choice of boundary conditions for quantum graphs resulting from con- 
sidering shrinking manifolds approching the quantum graph |10] ; no similar 
results are available for hybrid spaces, hence it is reasonable to take into ac- 
count a possibly large class of couplings. Our paper is based heavily of the 
the PhD-thesis of the second named author [16] where a class of the so-called 
disjoint boundary conditions was discussed. Later in [TO] another class of 
boundary conditions was suggested. We are treating here the most general 
boundary conditions. It comes out that this requires a certain extension of 
the notion of pseudoasymptotic expansion introduced in [18], but the most 
essential features of the coefficients and their interpretation are preserved. 

Acknowledgments. As already mentioned, the paper is based and extends 
the PhD-thesis of Svetlana Roganova. She would like to thank her advisor 
Jochen Briining for suggesting the problem and for many useful discussions, 
hints and comments. 

2 Definition of the Laplacian on a hybrid manifold 

2.1 Basic facts on self-adjoint extensions 

We recall here some basic facts on the theory of self-adjoint boundary con- 
ditions; for a more detailed presentation see e.g. [5] [6]. 

Definition 3. Let S be a closed densely defined symmetric linear operator 
acting on a Hilbert space H. Let Ti,r2 be two linear mappings from T>(S*) 
into a Hilbert space Q . The triple (G, T\, Y2) is called a boundary value space 
for S if the following three conditions are satisfied: 
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• for all x, y £ V(S*) there holds 

(x,S*y)-(S*x,y) = (F 1 x,F 2 



• for any u,v G Q there exists x £ V(S*) such that F\x = u and F 2 x = v, 

• kerTinkerTs = V(S). 

The existence of a boundary triple in the above sense is equivalent to 
the fact that the deficiency indices of S are equal, e.g. that S possesses 
self- adjoint extensions. Boundary triples deliver a useful machinery for de- 
scribing all self-adjoint extensions of S. 

In what follows in this subsection let (Q, Fi, F 2 ) be a boundary triple for 



Proposition 4. (1) Let A and B be bounded linear operators in Q satisfying 



then the restriction H A ' B of S* to the vectors f G V(S*) with ATif = BF 2 f 
is a self-adjoint extension of S. Moreover, each self-adjoint extension of S 
can be represented in such a form. 

(2) There exists a bijection between unitary operators U in Q and self- 
adjoint extensions H u fo S given byU <r+ H u = H A > B with A=l-U and 



Due to the above proposition the operator Hq which is the restriction 
of S* to kerTi is self-adjoint. For z ^ specHo the restriction of F\ onto 
ker(S'* — z) is a bijection; denote the inverse operator by 7(2). Clearly, 7(2) 
acts from Q to H, and we will refer to this operator as to 7 -field for the 
boundary triple, and the operator Q(z) := F 2 ^(z) will be called the Weyl 
function. 

Proposition 5 (Krein resolvent formula). For z ^ specH A ' B U specif 
there holds 



(H A ' B - z)' 1 = (H - z)' 1 - >y(z){BQ(z) - A) 'b^z)*. (2.1) 



2.2 Example 1: Laplacian on a segment 

In the Hilbert space Tt = L 2 [0,l}, I > 0, consider the operator S = —d 2 /dt 2 



with the domain V(S) = {/ € H 2 [0, 1} : /(0) = f(l) = /'(0) = /'(/))}• One 



can easily show that the adjoint operator S* is given by the same expression 
on the domain T>(S*) = H 2 [0, 1]. By direct calculation one can verify that 
as a boundary triple one can take 



S. 




B = i(l + U). 
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The distinguished self-adjoint extension Hq of S, i.e. the restriction of S* 
to kerTi is exactly the operator —d 2 /dx 2 with the Neumann boundary con- 
ditions. 

To construct the corresponding 7-field 7(2) and the Weyl function Q(z) 
it is useful to use the Green function Go(x,y; z), i.e. the integral kernel of 
the resolvent (Hq — z)^ 1 , z ^ spec£fo> which is given by 

cosh (\f^z(x — I)) cosh (^/— ~zy) 

— , x>y, 



, .^sinh (V—zl) 

cosh (y/—z(y — 0) cosn (V~^ x ) 

— , x<y, 



— zsinh (-y/— i/) 

where the branch of the square root ix fixed by ^VX > 0. In terms of Go 
one has 



l(z) (I) =H 1 G (;0;z)+t; 2 G (;l;z) 



and 

Q(z) 



G (0,0;z) G (l,0;z) 
G o (0,l;z) G (l,l;z) 



1 / COSh ly/^Z 1 



y/^zsinhly/^z V 1 coshZ^/-^, 
2.3 Example 2: Laplacian on a manifold 

Let M be a two-dimensional closed manifold and A = Am be the clo- 
sure of the positive Laplace-Beltrami operator defined initially on C°°{M). 
It is known that T>(A) C C°(M). Fix some (mutually distinct) points 
Qi ) • • • ) Qfj, £ M and denote by S the restriction of A to the functions vanish- 
ing at all these points (which makes sense, as all functions from the domain 
are continuous). The defficiency indices of S are (/x,/x). 

Let G(x, y; z) be the Green function of A, i.e. the integral kernel of the 
resolvent (A — z)^ 1 , then, for z ^ spec A, there holds 

V(S*) = {f = fo + ^jGCw) ■ fo e V(A), $ = (^-,^)eC}. 

The Green function can be presented as follows: 

G(x, y;z) = In 1 + F(x, y; z), 
2vr d{x,y) 

where the second term F is continuous in M x M. Using the above decom- 
positions one can show that there exist uniquely defined linear functionals 
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dj, bj on V(S*), j = 1, . . . , fj,, such that 



1 1 

= 2?r ln dfx y) +b iV> for x ~> 

and as a boundary triple for S one can take (C M ,ri,r2) with 

ri/ = M/)) J=w , r a / = fe(/)). =w 

Note that the distinguished extension of S corresponding to the boundary 
condition T%f = is just the initial operator A. By direct calculation one 
can show that the 7-field 7(2) and the Weyl function Q(z) for the above 
boundary triple are given by 



G(qj,q k ;z), j^k, 
F (qj,Qj; z ), j = k. 



More details to the above constructions can be found in [6]. We also mention 
the paper [12] treating a trace formula for the operators of the above type. 

Below we use the following asymptotic expransion, see [2J Formula (38)]: 
for all Af > 1, 

F(x, *,z) = ± (- In , 2 - 2 7 + £ r(n) ^' X) ) + 0(z-^). (2.3) 

\ n=l / 

Here a n {x,x) is the n-th local heat kernel coefficient and 7 is Euler's con- 
stant. 



2.4 Hamiltonian on a hybrid manifold 

Using the formalism of the previous section, we are going to construct 
Laplace type operators on hybrid manifolds and study some of their prop- 
erties. 

Consider a hybrid manifold constructed with m compact manifolds 
Mi, . . . , M m , with marked points qi s £ Mi, s = 1, . . . , and n segments 
Li, . . . ,L n , as described in the introduction. On each Mj consider also the 
usual Laplace-Beltrami operator A^, which is a self-adjoint operator on 
L 2 (Mi), with domain 2?(AmJ, the second Sobolev space of M%. Denote by 
Di restriction of Ajv^ to the domain 

P(A) = {/ € V(A Mi ) : f(q is ) = 0, 8=1,..., m}. (2-4) 

We parametrize the segments Lj by {xj S K : Xj G [0, Zj]}, where j = 
1, . . . , n and denote by D s - the closure in L 2 (Lj) of the operator —4^s defined 
on C^(Lj). 
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The operator S defined by 



s = Dx © • • • e D m © d{ © • • • e D s n (2.5) 

is a symmetric operator in L 2 {Mi) © • • • © L 2 (M m ) © L 2 (Lj) © • • • © L 2 (L n ) 
with deficiency indices (4n,4n). 

Definition 6. A Laplace type operator on a hybrid manifold is a self-adjoint 
extension of the operator S. 

We use the constructions of the previous section to describe all self- 
adjoint extentions of S. Namely, for each operator Dj take the boundary 
triple (O, Ti,Mji ^2,Mj) defined in subsection 12.31 and for each operator 
D 1 - take the boundary triple (C 2 , F\^. , ^2,Lj) as defined in subsection 12.21 
Then, obviously, their direct sum (C 2N , Fi, F2), where Fj = F j^@F j^, and 
T 3,M = Yj )Ml © ^j,M m , ^j,L = F jtLl © • • • © F j)Ln , j = 1,2, is a boundary 
triple for S. The corresponding 7-field 7(2) and the Weyl function are also 
the direct sums of the corresponding objects of each part. 

Below we use the enumeration of the points {qi S } ='■ {qj} in such a way 
that the point qj is the one attached to the initial point of if j = 2k — 1, 
k € N, or the the terminal vertex of if j = 2k, k £ N. 

Now we are going to restrict the class of admissible boundary conditions. 
It is natural to suppose that the boundary coniditions at each point of gluing 
are local, i.e. involves only the boundary values of the functions at this 
point. In view of the above enumeration of the points qj this means that we 
consider only the boundary conditions of the form AF± = BF2 where each 
of the 2N x 2N matrices A and B consists of four diagonal blocks 

A — ( ( ai '^ ( a i,i+N) \ B — ' ^>») (t>i,i+N) 



S a i+N,i) {di+N,i+N) J ' \(bi+N,i) Q>i+N,i+N] 

It is useful to introduce the 2x2 matrices 

A - — I ai >* a i,i+N \ £> . _ I ^i,i+N 



clearly, these matrices should satisfy the same conditions AjB* = AB* and 
det(AjA* + BjB*) / for each j. It is also clear, by proposition 2J that 
each of the boundary conditions can be rewritten using unitary matrices of 
the same structure 

Aj = l-<f>j, B j = i(l + ^ j ) 

with 

\<Pi+N,i Wi+N,i+N J 
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and 

\{<Pi+N,i) {(Pi+N,i+N) J 

This paramatrization will be referred to as canonical. Note that diagonal 
matrices <J?j correspond to the absence of coupling at qf the boundary 
conditions split into one boundary condition on the manifold and another 
one on the segment attached; we will call such boundary conditions reducible. 
Below we consider non-reducible boundary conditions only, hence assuming 
that all the matrices $j in the above constructions are non-diagonal. 



3 The resolvent expansion 

Warning. From now on, due to technical reasons and in order to be com- 
patible with other works, we will use —z 2 as our spectral parameter, and not 
z. However, we will still use e.g. the notation Rq(z) to denote (So + -2 2 ) -1 , 
and not (So — z)^ 1 as in the previous section. Note for further use that in 
this notation the Green function on manifolds or segments obeys 

G(x, u, z)G(u, y, z) du = -— G' z (x, y, z) 

M lz 

(which is just a consequence of the Hilbert resolvent identity). 
3.1 Computation of Tr R 2 

Let us start with the computation of Tr R 2 using Theorem Setting 

T{z)= 1 (z)[BQ{z)-A]- 1 B{ 1 {z))\ 

we obtain 

Tr R 2 {z) = TtRq(z) - Tr(R (z)T(z)) - Tr(T(z)R (z)) + TrT 2 (z). 

We recall that the resolvent Rq is the direct sum of the resolvents of 
the ordinary Laplacians on the manifolds Mi (denoted by i?o,Mi) and the 
resolvents for the Neumann Laplacians on the segments Lj (denoted by 
Ro t L ) forming the hybrid manifold. For the manifolds one has 



Tr Rl M .(z)= / G Mi (x,u,z)G Mi (u,x,z)dxdu 

J Mi J Mi 

-{G Mi )' z {x,x,z)dx, 



'Mi J Mi 

1 



Mi 2z 

and for the segments one uses the same formula with the explicit expression 
(12.21); we will use it later. 
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The linear operator T(z) can be also rewritten as an integral operator 
as follows. Denote the entries of the 2N x 2N matrix [BQ(z) — A]~ 1 B by 
Cij(z). Using G(x,y,z) = G(y,x,z) we get, for any / G L 2 : 

T(z)f(y) = j(z)[BQ(z) - A^B^Tfiy) = 
^2 / c ij (z)G(y,q i ,z)f(u)G(q j ,u,z)du. 

i,3 J 

The integral here is in fact the sum of the integrals over manifolds and 
segments. Using this integral representation of the operator A, we calculate 
the remaining terms of Tri? 2 as follows. First, for any / G L 2 : 

R (z)T(z)f(x) = Y^Cij( z ) J J G(x,y,z)G(y,qi,z)f(u)G(qj,u,z)dydu. 



1,3 



Hence, the operator Rq(z)T(z) has the integral kernel 

K(x,t,z) = Y^Cij(z) J G(x,y,z)G(y,qi,z)G(qj,t,z)dy, 



1,3 

therefore 



TrR (z)T(z)=Y^Cij(z) J J G(x,y, z)G(y,qi, z)G(qj , x, z) dx dy, 



1,3 

or, 



TtR (z)T(z) = ^Y^Cij(z) J (^j G(x,y,z)G(y,qi,z)dyj G(qj,x,z)dx 
+ \Yl Ci ^ j (/ G{x,y,z)G(qj,x,z)dxj G(y,qi,z)dy 
= \^ c ij( z )( j -7^ G ' z {x,q i ,z)G{q j ,x,z)dx+ 

i,3 

J -^G' z (qj,y, z)G(y, q iy z) dy^j 
= ~^Yl Ci ^i z ) J (G(qj,x,z)G z (x, qi ,z) + G' z (q j ,x,z)G(x,q i ,z)) dx 

i,3 

= ~\^Yl Ci 3^ J (G(qj,x,z)G(x,qi,z))' z dx 

i,3 



g - z , ' ' ' ' \ - G . [(ij.q,. z 



Cij[Z) 



1,3 



G' z ( qj , qi ,z) G'^q^zy 
8z 3 8z 2 
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Let us now find Tr T 2 (z). We will use the same method and represent 
the operator T 2 (z) as an integral operator: 

T 2 ( z )f = X^-W / G ( x 'Hi z ) G (ij^y^ z ) T ( z )f(y) du 

1 z^f^^ 



i,j,k,l 

G{y,q k ,z)G(q l ,u,z)f(u)du) dy 



= ^ Cij(z)cki(z) j G(x,q i ,z)G(y,q k ,z)G(q j ,y,z)G(qi,u,z)f(u)dudy. 
i,j,k,l 

Now we find: 

i,j,k,l 

= ^2 ^ c i j(z)c k i(z)G' z (qi,q i ,z)G' z (q j ,q k ,z). 

i,j,k,l 

Summarizing the results obtained, we state 

Theorem 7. Consider a hybrid manifold H , consisting of the manifolds Mi 
and the segments Lj, and the Laplace operator corresponding to the boundary 
conditions ATi = BT2 on H. Then the following formula for the trace of 
the square of the resolvent of this operator holds: 



Tr R 2 (z) = — I —G' z (x,x,z)dx 
Jh % z 



2N 



G' z (qj,qi,z) G" z (qj,qi,z) 



2N 



8z 3 8z 2 
+ Z~2 Yl c iA z ) c kl(z)G' z {ququz)G' z (q^q k ,z). 



Az 2 

i,j,k,l=l 



Here G(x, y, z) is the Green function of Sq on the hybrid manifold and the 



entries of the matrix [BQ(z) — A] l B are denoted by (Hj(z)i<i j<^2N- 



3.2 The asymptotic expansion 

The result of Theorem [7] is given in terms of the Green functions for the 
Laplacians on the smooth parts of the hybrid manifold. But using it in 
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applications is practically impossible for two reasons: first of all, inverting 
the 2N x 2N matrix BQ(z) — A can be difficult, and secondly, an explicit 
computation of the Green function is almost never possible. Nevertheless, 
there are ways to get some simplifications if we restrict attention to large z. 
As we will see below, the special structure of the matrix (BQ(z) — A) allows 
us to find its inverse asymptotically as z — > oo. Likewise, one can use the 
representation of Rq as an asymptotic series in powers of z (analog of heat 
kernel expansion). Thus we may attempt to generalize (|1.3p to the case of 
hybrid manifolds. 

For further calculations the following lemma is crucial. 

Lemma 8. Let M be a compact two-dimensional manifold and G(x, y; z) 
be the integral kernel of (A + z 2 )^ 1 . For any x,y G M, x 7^ y, and any 
k = 0,1,2, ... , there exists C k = Ck jX ,y > such that for z — » +00 one has 

d k G(x,y,z) =0(e -C kZ) 
dz k 

Proof. Let P(x, y; t) be the heat kernel of A, i.e. the integral kernel for 
the semigroup e~* A . It is known [B, Corollary 6] that there are positive 
constants a and 5 such that 



0<P(*, y; t)<%xp' d ^ V? 



Clearly, 



hence 



t *V 4(l + <J)t 
/•+00 2 

G(x,y,z) = / e~ z 'P(x,y;i)cft, 
J 



d k G(x, y, z) 



p k {z,t)e z t P(x,y;t)dt, 



dz k 

where pk is a certain polynomial. One can estimate 

d k G(x, y,z) ^ ^ f + °° q k (z, t) ( d(x, yf \ __ z * t 



<a v exp - -A ,a/ ) e - z l dt, 

~ 1 t L \ 4(1 + 5)t J 



dz k 

where q k is another polynomial. Now let us note the identity 

r+00 

/ t^e-^e-P'dt = 2(b/p) v / 2 K v (2^/bp), 
Jo 

where K u is the modified Bessel function, see [151 Supplement 4.3], which 
gives, with some polynomials qj, 

d k G(x,y,z) A fd(x,y) 2 y/2 /zd(x,y) 

<2a^ qj (z)(——) K 



dz k 



4(1 + 6) ) " J \VT+S 



Now it is sufficient to recall the asymptotic behavior K u (s) = 0(e s ) for 
s — ► +00. □ 
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Thus, if the points x and y do not coincide, the Green function G(x, y, z), 
as well as its derivatives G' z (x,y,z) and G zz (x,y, z) decay exponentially as 
z tends to infinity. 

In the statement of Theorem [7] the terms G(qi,qj, z) appear, where the 
distance between the points qj and qj is fixed by the configuration of our 
hybrid space. Hence all terms of such type for non-coinciding points qi 7^ qj 
are exponentially small as z —* 00, so we need to take into account only the 
diagonal terms G(qi, qi, z). Using this observation and performing calcula- 
tions similar to the proof of Theorem but neglecting terms of exponential 
small order of z we prove that for some positive constant c, as z —* 00 

TiR 2 (z) = — I —G' z (x,x,z)dx 

, 9 V 1 ^ ( G 'z(QuQu z ) G'!, z (qi,qi,z) \ 

i=l x ' 

. 2N 

+ 4^2 ^2 cij^cjii^G'ziqiiqii^G'ziqjiQj' 2 ) + o{e~ cz ). 

This result can be rewritten in a compact matrix-form, as we are going to 
see now. Let us denote by G' z the diagonal matrix with entries G' z (qi, qi,z). 
Denote also by G" zz the diagonal matrix with entries G zz (qi,qi, z). This is a 
slight abuse of notation, since G' z and G zz are also well defined operators; 
it will be clear from the context what is meant by this notation. Then it is 
easy to see that Q can be formulated as follows. 

Theorem 9. We have the following asymptotic relation as z — > 00: 
Tr R\z) = -J h dx-^Tr {G" zz [BQ - A]^B) 

+i ^3 Tr (G' z [BQ - A^B) + -L Tr (G' z [BQ - A^B) 2 + 0(e~ cz ). 

Now we need only to compute (BQ(z) — A)~ l B, and fortunately, this is 
possible in this approximation. 

As we have already mentioned, the matrix Q(z) has a direct product 
representation, Q(z) = Q\ @ . . . Q m © G where Qi is the Q-matrix for the 
i-th manifold Mj in the hybrid manifold (see subsection 12. 3p . for example, 

/ F(qi,qi,z) ... G(qi,q^,z) 

Qi = 

and G = (Gki) is an TVx A^-matrix consisting of the Green functions (denoted 
by G s ) for the Neumann Laplacian on the segments: 

G s (qk, qi, z) if qk and qi belong to the same segment, 
otherwise. 



G 



kl 
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Nevertheless, the matrix BQ(z) — A is still too complicated to find its 
inverse explicitly. But letting z — > oo and applying Lemma [8] the situation 
improves radically. More precisely, all non-diagonal terms of Q(z) are either 
zero or of type G(qi,qj, z), i ^ j, and decay exponentially as z tends to 
infinity. Therefore, we can write Q(z) = Qd(z) + Q(z), where Q(z) is ex- 
ponentially small as z goes to infinity, and the leading order approximation 
Qd(z) to Q(z) is given by 

Qd{z) = diag (^F(q 1 ,q 1 ,z), . . . , F(q N ,q N , z),G s (q 1 ,q 1 , z), . . . ,G s (q N ,q N , z)j. 

Then (B(Q{z) — A)^B to first order, i.e. (BQd(z) — A) , consists of four 
diagonal N x N blocks. The inverse matrix can be found using the Frobenius 
formula for a block matrix consisting of matrices A, B,C, D: 

A By 1 _ ( (A- BD^C)- 1 A^BiCA^B - D)~ l 
CD) ~ \{CA- l B - D)- l CA- 1 (D-CA^B)- 1 



Hence one can obtain the block representation 

[BQd(z) - A)^B - 
where U, W, V, X are diagonal N x N matrices with 



x- l w x- x v 



Xi^i = Xi = {bi^bi-\-N,i+N — bi,i+Nbi+N,i)FiGi 

+ {ai+N,ibi,i+N — ai,ibi+N,i+N)Gi + (aij + Nbi + N,i — ai+N,i+Nbi,i)Fi 

+ {ai,i(ii + N,i+N — ai,i+Ndi+N,i), (3.1) 

and 

Ui t i = Ui = (bi^bi+N^+N — bi,i+Nbi+N,i)Gi + (cii : i + Nbi+N,i — a,i + N,i+Nbi,i), 

Wi^i = Wi = di^ + Nbi+N,i+N — ai + N,i+Nbi,i+N, 

V%,i =Vi= (b iti bi + N,i+N — bi,i+Nbi+N,i)Fi + (di + N,ibi,i+N — a*,ibi+N,i+N)- 

(3-2) 

Note that if we change the parameters A — ► LA, B — > LB, with a non- 
degenrate matrix L of the same structure as A or B, then all of the above 
expressions acquire the factor det Lj . 

It is important to mention that the matrix Xj is invertible for large z. 
To see this, let us note first that all four summands in (|3.ip have different 
asymptotics for large z: by (|2.3p one has F{ ~ Inz 2 , and by direct calcula- 
tion, 

G i (z) = l/z + 0(e- cz ), 

hence it is sufficient to show that at least one of them is non-zero. To 
show this, it is useful to assume that the matrices A and B are taken in 
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-1 



the canonical form, see subsection 12.41 In this case, assuming that all four 
summands in (|3,ip vanish one arrives for each % to the incompatible system 

1 + <Pi,i + <Pi+N,i+N + det $i = 0, 1 - tpij - <fi + N,i+N + det = 0, (3.3) 
1 - <Pi,i + <fi+N,i+N ~ det $j = 0, 1 + - (fi+N,i+N ~ det $j = 0. (3.4) 

Lemma 10. With the notation above, we have 

[BQ(z) - A}- 1 = [BQ d (z) - A]- 1 + 0{e~ cz ), 

for some positive constant c. 

Proof. First, we write 

[BQ{z) - A}- 1 = [B(Q d (z) + Q(z)) - A 

= [BQ d (z) - A]- 1 ^ + BQ(z)[BQ d (z) - A]' 1 ]' 1 . 

As can be easily seen from the expressions for Uu, Wu and Va, one has 
[BQd(z) — A]~ 1 B = 0(z 2 ). As Q(z) decays exponentially fast for large z, 
the product [BQ d (z) — A]~ 1 Q(z) will then also be exponentially small. By 
using the Neumann series to compute [1 + [BQ d (z) — A]~ 1 Q(z)]~ 1 , we see 
that this term is of the form 1 + 0(e~ cz ), and this immediately implies the 
conclusion of the Lemma. 

□ 

Lemma [TU] allows us to reformulate theorem as follows: 
Theorem 11. For z — > +cx> there holds, up to 0(e~ cz ) terms, 

TlRHz) = _[ ote,,,,) A _ i y-irpi+gry. 

J H 2z Az ^ Xi 

1 Fp^Gfr 1 {FlU l f + 2F' l G[\W i \ 2 + {G' i V l f 
+ 4z3Z^ Xt + 4z2Z^ X f ■ {6 ' b) 

i i 1 

This theorem describes the large spectral parameter asymptotic of the 
trace of the second power of the resolvent for a Laplace operator through the 
Green functions of the manifolds and segments forming the hybrid manifold. 
In the next section, we will write this expansion in terms of heat kernel 
coefficients for the manifolds forming the hybrid manifold. 

3.3 Pseudoasymptotic resolvent expansion in terms of heat 
kernel coefficients 

To compute the above expression for the trace in terms of heat kernel coef- 
ficients, we need to consider the term 



G' z (x,x,z) 
2z 



15 



In fact, this term consists of two types of integrals: over the manifolds and 
over the segments. As we already know (see equation (jl.2p ). for manifolds 
the following expansion holds as z — > oo: for all q ^ 

" J ^tr^ dx = Tr ^ = t aj §^ + o(z-^), 

where Rq is the resolvent of the Laplace operator on this manifold and the 
constants are the heat kernel coefficients for it. The integral over segments 
can be easily computed asymptotically using the explicit expressions (j2.2j) : 



2z 4z 3 2z 4 

Definition 12. Assume that there exist rational functions Cf~, k = 1, 2, ... , 
such that for any q > and a £ (0, 1) a function f = f(z) admits the 
following representation for z — ► +oo: 



k=0 v 7 



(3.6) 



An expansion of the form (|3.6p will be called pseudoasymptotic, and in this 
case we write 

fc=0 

Lemma 13. If a function f admits a pseudoasymptotic expansion, then this 
expansion is unique. 

Proof. First, we notice that if P and Q are two polynomials and a > 0, 
then, for large z, 

\Q(lnz 2 ) 

Suppose that there are two pseudoasymptotic expansions of /, i.e. for any 
q > and any a > 

m--±^ + o{^)--±^l + o{^). ,,r, 

fc=0 V 7 fc=0 v 7 

For q = we obtain 

c (lnz 2 ) -d (lnz 2 ) = 0(z~ Q ) 

which implies Co = do- By an immediate induction on k, we show similarly 
that we will always obtain c^lnz 2 ) — G?fc(lnz 2 ) = 0(z~ a ) and consequently 
the uniqueness of the expansion. □ 
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Remark 14. Let us notice that if a function / possesses a pseudoasymptotic 
expansion, and if each coefficient is a rational function with nonpositive 
degree, then each Ck can be uniquely expanded as a In z 2 -asymptotic series: 
for all K > 0, 

for some constants c^i- 

Theorem 15. The trace Tr R(z) 2 obtained in Theorem [771 admits a pseu- 
doasymptotic expansion of the form 



47TZ 2 4z 3 ^ z fc 

fc>4 

EjX(Mj) A A„ , 1 



where 04(A) = — —— h — H — - — h of — I , A — ► +00 and 

o 4 A V A / 

N o = #{j G {1, • • • ,N} : a iti+N b i+ N,i - a i+ N,i+Nk,i °}- 

Moreover, if Nq = N , then c& = O(l) as z — > +00 for all k. 



Proof. Consider first more carefully the expression (|3.ip for X{. Let us show 
first that at least one of the last two terms on the right-hand side, a.jF{ with 

otj = ai,i+N0i+N,i — a<i+N,i+Nbi,i and (3j := (ii,i(H+N,i+N — di,i+Na>i+N,i = 
det A,- is non- vanishing. Taking the matrices Aj and Bj in the canonical 
form An = 1 — Uj, Bj = i(l + £/,•), and assuming that the two above terms 
vanish one arrives at the system 

1 + (p iti - (p i+N>i+N - det $j = 0, 1 - </3»,i - (p i+N)i+N + det $j = 0, 

and one immediately obtains ifi+N,i+N = L which means that the matrix 
<£j is diagonal. As this case is excluded by assumption, we have the result. 
Hence on has always 1/X, = 0(1) for large z. 

Consider the second sum on the right-hand side of (j3.5|) : the other terms 
are considered in the same way. Each summand is of the form 

_ Ff'Ui + G'jVi 

E > = Xi 

Assume first that ai ; i + Nbi+N,i — ai+N,i+Nbi,i 7^ 0. In this case, substitut- 
ing the expansion (|2.3|) for F{ into the expression in Theorem [TT] one arrives 
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at an expansion of the form 



E 



^5 I - 2^fc=0 z fc ^K+l 



K 

1 ' IS 

lnz 2 d' fe „, 1 



zln(zVo) z fc ln(z 2 e d o) 

For large z the last term in the denominator can be expanded as 

-i 



zln(z 2 e<>') ^z fc ln(>Vo) V ^ +1 

1 A ^ln Sl (z 2 ) 
= 1 + > — — - — h O 



2 L+1 



i ln S2 (z 2 e d o)^ \z L+1 

for some constant and non-negative powers si,S2- On can see also that 
si ^ S2- Multiplying both series in powers of z in the expression for the first 
fraction we obtain an expansion of type 

m = A« / i 
^ln^Vo'y U L+1 / ^ zl 

with coefficients 

Q;(ln z 2 ) 

which are rational in lnz 2 . Note that in this first case we always have 
degP; < degQi. 

The case a,i t i + ]ybi+N,i — CLi+N,i+Nbi,i = can be treated in a similar way, 
one deals with an expansion of the form 



E 



z 5 + 2^k=0 z k + u \ K+l 



3 ' In -y 2 r]" 

1 + ^ + EE*4+0(jAt) 



and develops the denominator directly using the geometric series. 

Furthermore, examining the asymptotic behavior of all the functions, one 
can easily conclude that the last three sums on the right-hand sum of (|3.5[) 
are of order 0(z~ 4 ), and the only terms that contribute to the coefficient C4 
come from the summands 

F'.'TJ- F'TJ- 

Kj := _4±Xl + h( 1 



Az 2 Xi Az 3 Xi 
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Note that if ai,j+jv^i+jv,i — a-i+N,i+Nbi,i = 0, then U{ = 0{z~ l ) and Kj does 
not contribute to C4. Otherwise, by considering the leading terms in the 
numerators and the denominators, 

K 

In z 2 
Summing over all j we obtain the result. 



□ 



Remark 16. It is useful to emphasize that all the "coefficients" c& are non- 
constant, i.e., we have always some logarithmic terms. Indeed, it is clear 
that the logarithmic terms can only disappear if the coefficients of the terms 
containing F{ in the expressions (|3.1|) and (|3.2|) vanish. 

Taking the matrices A and B in the canonical form (subsection I2.4p . this 
condition would mean that 



1 + <Pi,i + tpi+N,i+N + det $j = 0, 1 + - ip i+Nt i +N - det $j = 0. 

This implies directly tpa = — 1, and the matrix $j must be hence diagonal, 
which contradicts to the non-reducibility of the boundary conditions. 

Remark 17. The expansion obtained allows one to address some questions 
of the inverse spectral theory. It is clear that the knowledge of the coef- 
ficients in this expansion gives the summary volume of the manifolds and 
the summary length of the one-dimensional parts. Considering the fourth 
term one can obtain some estimates about the Euler characteristics of the 
manifolds and the number of the segments. Note that for generic boundary 
conditions one has Nq = N hence generically one can recover the sum of the 
Euler characteristics of the manifolds and the number of the segments. We 
will adress such questions in greater details in the next section. 

Let us first note that the formula for Tr R 2 obtained in Theorem [15] 
depends on the heat kernel coefficients for the smooth parts of the hybrid 
manifold, and it is well known that heat kernel coefficients are recursive. 
Some recursivity of the coefficients in the expansion of Tr R 2 follows as well. 
But as a series inversion is required, we cannot obtain this recursive formula 
in an explicit way. Nevertheless, it is possible to find the formula for some 
terms in this expansion for some special class of boundary conditions. 



4 Special case of boundary conditions 

As can be seen from the proof of theorem [T5l the boundary conditions enter 
the pseudoasymptotic expansion in a rather involved way. In this section 
we treat in greater details the case when the matrix B is invertible. In this 
case one can normalize by B = 1, and A must be then self-adjoint with the 
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same block structre as previously. To distinguish from the previous case on 
the notation level, we denote the matrix elements of A by A, 



A = ( (\i+n) 
\(\,i+N) (K+n,i+n) 

and the diagonal entries are real. We assume additionally that the entries 
Ai+jv^+jv are non-zero for all i G {1, ... , N}. 

Note that the expressions (|3.ip and (|3.2p simplify, 

■^Q, = (Ft ~ \,i)(Gi — \+N,i+N) — |Aj.i+Ar| 2 , 

Ui = Gi — Aj+7v,j+Af) W{ = A^j+AT, Vi = Fi — 

Theorem [PS] is, of course, applicable, and one has additionally Nq = N. As 
noted in the proof of Theorem 1151 we also have c(A) = O(l) for large A. 

To simplify notation in this section, we will use the following convention. 
If the function / = f(zAnz 2 ) has a z-pseudoasymptotic expansion with 
coefficients c q (lnz 2 ), we noticed before the statement of Theorem 1151 that 
each Cq(lnz 2 ) can be expanded as 

c ^ ln z2 ) = E + ° [ (ln z 2 )K +i 

for some constants c q k- We will then write 



/- E 



Cqk 



z gn n z 2\k- 
q,k=0 V ' 

4.1 First terms in the resolvent expansion. 

Lemma 18. Let n > 4. The coefficient c n appearing in Theorem{T^ has an 

expansion in 1/logz 2 , and the first term in this expansion is of the following 
form: 

Zixm , n n _ 4 . 

6 + 4 ' " ' 

^ 2k -I 

> n = 2k + l,k>l; 

i=l ^ A i+N,i+N 

Eak-;T(k + 1) v-^ 2k 
47r + E AX 2k~2 . n = 2k + 2,k>l. 

i ^ i=l i+N,i+N 

Proof. The terms in C4 and the first term in c n for even n arise from the 
expansions of the second power of the resolvents for all manifolds and seg- 
ments. The other terms require some calculation. Let us look once more at 
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the expression in Theorem 1111 and analyze it more carefully. As we can see 
the expansion of 



CFi-Ai)(Gi-A 2 )- |A 3 | 2 ) 

always has Inz 2 in the denominator, i.e. for some non-zero constants k'^ m 
one has the z-pseudoasymptotic expansion 

J_ 47r Km 

A ~ \ l+N) i + N^z 2 + n ^ o z«(lnz2)m+i> 

where we denote (Fi — Aj j j)(- — Aj+jv.i+iv) — |Ai,i+jv| 2 ) by Dj. In the corre- 
sponding numerators we have terms of type F", F[ and F{. Their expansions 
as z — > oo are: for all p ^ 1 

Ft(x, x,z) = ± (-2 7 - In + ^g^ ] + 0(z~^), 

\ n=l " / 

*?(», «,,) = ^ - ± 2nTi ^' X) + 0(^- 2, - 3) ). 

n=l 

„, J_ A 2n(2w + l)r(w)a tli (s,3;) 2p _ 4 

^ aX,X,Zj ~ 2vrz 2 + ^ 4vrz 2 ™+ 2 + 1 J ' 

n=l 

It is clear now that the first terms in the expansion of c n (i.e. terms without 
lnz 2 ) can appear only in terms containing Fj. All other terms will contain 
this logarithm in denominator. Moreover, only the part ~ l ^ JT z in Fi matters. 

The same arguments are valid for the last summand whose denominator 
contains (lnz 2 ) 2 , and consequently we should take only the term (Fi — A, i) 2 
into account. Finally, we try to find the terms containing only powers of z 
in 

— 2Fi — Fi (Fi — Xi^i) 2 



4z 5 Di 4,z e (Di) 2 ' 

First of all let us treat the denominator more carefully 

1 4tt f i 1 { y K 



A \ i+ N,i+N In z 2 y z\ i+N)i+N ^ z n In z 2 ^ 

a I OO OO j ft 



2\m+l I ' 



X i+NMN lnz 2 \ ^ (z\ i+N4+N ) n n ^ =(j z n (lnz 

with some coefficients k' n ,k'^, whose explicit expressions are not important 
now. 
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Then we see that in . f J the terms containing only powers of z (without 
logarithm) are 

-3 4tt ( ^ 1 \ -lnz 2 



4z 5 X i+N , i+N lnz 2 \ ^ (zX i+N>i+N ) n 



An 



3 A 1 



4 Z-^ y n+5\n+l 
n=0 Z A i+N,i+N 

Performing the same reasoning we find that 



(4tt) 2 ^ n + 1 ^ K 



27; (A i+JV , i+Jv ln^) 2 ^ ^ Z «A» W ^ ^(ln^)m+i j ■ 

The only term which cancels this logarithm in (Fi — Aj^) 2 is ^j^p - - And 
the contribution to the "pure" polynomial part in z is 

_L_ (4tt) 2 y> n + 1 (lnz 2 ) 2 _l^ n + 1 

4.6 • (A^.^lnz 2 ) 2 2- " (4vr) 2 " 4 2- z^X^. +N 

Finally, we find (not taking into account those which arise from the 
expansion of Rq) 



3 \ - 1 1 \ - n+1 1 x - n + 3 

I 2-~s n +5\n+l ' 4 2-~s n +6\n+2 ~~ 7 7J^7Tn+l ' 

* n=0 Z A i+JV,i+AT * n =0 Z A i+JV,i+AT * n=0 Z A i+N,i+N 

Collecting the terms with the same power of z we prove the lemma. □ 

The next result requires more complicated calculations of the same na- 
ture, so we just state the main steps of its proof. 

Lemma 19. Fix some integer n > 4 and consider the coefficient c n which 
appears in Theorem [731 The second term of the expansion of c n has the 
following form: 



nyn -^-^ +m n = 2k+lk>1 . 
i=l A i+N,i+N 111 z 



^ ir(n - A)(n - 2)\Xj, i+N \ 2 ^ a H (l + 1)! 

E 7^=3 7— 2 + E i QZ 2 > n = 2/ + 4,/^0. 

Proof. The idea is the same as before: one can obtain the terms of this form 
only from specific terms in the expansion. Due to the special structure of 
the denominator we have terms of type and terms g ^V^ an d the fact 
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that in the numerator only the terms with Fi contain a logarithm, allows us 
to restrict our calculation to some specific terms. 

We need as before the expansion of tt but up to the second power of In z 2 
and the expansion of-^s up to the third power. To simplify the expressions 
we denote the following quantity by Wj 



i+N,i+N J Z \+N,i+N 
oo 



1 \ (n — l)!«r. 



+ — E 

V zXi+NA+N J ^ Z 



i+NA+N 



7 2n 



n=l 



and state that up to the second power of In z 



2 



1 47T ^ ( I 



A ^i+N,i+N In Z 2 ^ \ zAj + Ar,i+AT 

4ttw;- ^ / i x n - x 

n 

n=l 



Xi+N,i+N hi 2 Z 2 ~ V zXi+NA+N 



and up to the third power 

1 16vr 2 ~ n 



E 



A 2 A 2 +7ViJ+iV In 2 z 2 ^ (^^i+Jv)"- 1 
A-K 2 Wi ^ n ( n _ i) 

Then the terms we are looking for appear in 

Di y Az 2 4z 5 Az 3 

and in 



1 ({Fi-X 



2 

i.i J 



D 2 V 4z 6 

Expanding and summing these expressions we find that the terms con- 
taining the first power of logarithm in the denominator can be arranged in 
two sums as in the statement of the Lemma. □ 

In the same way, but using much more complicated calculations, we can 
get the following result whose proof we omit (see |16] for details). 
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Lemma 20. In Theorem [75| the terms with In 2 z 2 in the denominator of 
the c n are (we do not separate the different powers of z in order to not 
complicate the expression): 

1 A/l 2 7 4vrA 2i ^ (2n + l)(n-l)!a ni , ,„ . 

^2.1^-^-^ + 2. ^rz + ( 2 ^ 

j=l \ n=l 

2vrA| i ^ n(n + 2) ^ 2(n + l)n\a ni 



y 2n+4 



n=l n=l 

27r 2 A3j ^ n ( n 2 _ i) ~ ( n + l) n \a ni 47rA|- ^ n + 1 

±l n=l n=l n=0 



n=l fc=l 

oo 



(to- l)!a ni ^ (k + l)/c!a^ 7rA| / v Mn - l)!a n i ^ fc(fc + 2) 

+ 2. y 2n 2. z 2fc+4 + Z 4 A 12. z 2n 2.( zA U 

x n=l fc=l 



E2(2n + l)n!a n j 1 2n!a n j 1 + 2k 

i 2 ^ 2w ( 2 A H ) fc + ^ "i 2 ™" ^ (zAii)* 

n=l fe=0 V ; n=l k=0 V U ' 

where \u = \i+N,i+N, = K,i and X^i = \\i t i+N\. 



4.2 Inverse spectral data 

Being restricted to the class of the boundary conditions introduced in the 
beginning of the section, we consider now the inverse spectral problem, i.e. 
the question "Which information about the initial system can one obtain 
using the expansion of the second power of the resolvent?" . As we will see, 
it is possible to find some geometric characteristics of the system and some 
information about the operator on it. 

Theorem 21. Consider the expansion of the trace of the square of the re- 
solvent as in the Theorem \15[ The knowledge of TrR 2 determines: 

• whether this manifold is hybrid or smooth; 

• the sum of the volumes of all manifolds taking part in the construction; 

• the sum of the Euler characteristics of all manifolds; 

• the number of segments used in this hybrid manifold; 

• the sum of the lengths of these segments; 

• the Euler characteristic of the hybrid manifold. 

Proof. The presence of In z 2 -type terms is a criteria of singularity. If in the 
expansion of TrR 2 there are no such terms, this means that the considered 
manifold is a "normal" manifold without any singular points. Indeed, the 
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log-terms appear only from the singularities of the Green function (and not 
its derivatives) at the points of gluing. 

The coefficient of z~ 2 is equal to ^volMj/(4-7r) and provides us with 
the sum of the volumes of all manifolds. 

The coefficient of z~ 3 is equal to ^£j/4 and provides us with the sum 
of the lengths of all segments. 

Considering the term of type l/(z lnz 2 ) we find the number N/2 of all 
segments. 

The coefficient of z~ 4 is equal to 

Ex(Mj) N 
6 4 

and gives as the sum of the Euler characteristics of all manifolds since we 
already know the number N. 

Finally, using a Mayer- Vietoris type argument, see [16] for details, it is 
not difficult to show that the Euler characteristic of the hybrid manifold is 
x{Mi) — N, which is easy to find now. □ 

The results obtained in this work do not let us find the volume of each 
manifold and the length of each segment separately. But it should be possible 
to find these quantities with the use of scattering theory [14] . 
The resolvent expansion provides us also with some information about how 
we glue the segments to the manifold in the hybrid space. Namely, we can 
obtain some information about the matrix of boundary conditions which 
defines the Laplace operator on the hybrid space. 

Theorem 22. Consider the z-pseudoasymptotic expansion of the trace of 
the square of the resolvent expansion of type [731 If we assume that we know 
the heat kernel coefficients for all manifolds composing the hybrid manifold, 
and that the coefficients \i + N,i+N are mutually distinct and nonzero, we can 
find the diagonal elements of the matrix of boundary conditions A and the 
absolute values of its non-diagonal elements up to permutation. 

Proof. As it was shown in Lemma [18] the first terms in asymptotic expansion 
of c n are of the following form: 

Em 1 x(m 1 ) n 

6 + T' n = 4; 

TV 



2k - 1 

-^Tfi , n = 2k + l, k> 1; 



i=l ^ A i+N,i+N 



EakiT(k + 1) x— > 2k 



Air ^-^ A \ 

M x i=l i+N,i+N 
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Consider Xi = t , i = 1, • ■ ■ , N. Taking the first N coefficients of 

powers of z starting from z 5 and denoting them by cZ ri _|_4 , n = 1, • • • , N we 
obtain iV equations which have the foilowing structure 

d n +4 = X i 

i 

To find the soiution of such a system we find the following symmetric poly- 
nomials in terms of c n using a well-known recursive procedure: 

5*1 = x\ H h xn, 

S 2 = X\X 2 + X1X3 -\ h X2X3 + . . . , 

S3 = XiX 2 X 3 + XlX 2 X 4 -\ 1- X2X 3 X4 + . . . , 

Sn = x\X2 ■ ■ ■ xn- 

Due to the Viete Theorem we state now that x±, . . .xn are the roots of the 
equation 

x N - S 1 x N - 1 + S 2 x N - 2 + ... + S N = 

and can be found up to permutation because we assume that they are mu- 
tually distinct. 

If we substitute the obtained values of Aj + jv,i+iV in the expression ob- 
tained in Lemma \W\ we obtain a system of linear equations for lA^j+Arj and 
can find them. 

As soon as we have found these elements, the result of Lemma [20] gives 
us a system of linear equations for A^j and we can also find them. The 
Theorem is proved. □ 
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